Electron transport in ZnO thin films 
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Epitaxial, n-type ZnO films grown by a laser molecular-beam epitaxy method were investigated by the 
temperature-dependent Hall-effect technique. The 300-K carrier concentration and mobility were about 
n s ~ 10 16 cm -3 and 440 cm 2 /Vs, respectively. Transport characteristics are calculated by solving the Boltz- 
mann transport equation using a variational method. Mobility limit of 430 cm 2 /Vs was calculated at 300 K. 
The temperature dependence of the mobility for an undoped film is calculated and agrees favorably well with 
experimental data if physical parameters are chosen so as to approach those. In the experimental 'mobility 
versus concentration' curve, unusual phenomenon was observed, i.e., mobilities at n s ~ 5x 10 18 cm -3 are 



significantly smaller than those at higher densities above ■ 
electron-plasmon interaction. 

PACS numbers: 78.55.Et: 81.15.Fg: 71.35.Cc: 72.15.-v 



10 cm . It is qualitatively explained in terms of 
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ZnO is direct band gap semiconductor suitable for the de- 
velopment of efficient blue and UV optoelectronic devices. 
Attempts have extensively been made to incorporate shallow 
acceptor impurities in this material. We have recently devel- 
oped a reliable and reproducible method to grow p-type ZnO 
layers, leading to the observation of electro-luminescence 
from a homoepitaxial p-n junction structure Q]j. The exploit 
for discovering the p-type conduction has noticed us the im- 
portance of quality improvement of undoped epilayers as a 
starting point of reliable doping Q], In this work, the elec- 
trical properties of epilayers are investigated since mobility is 
considered to be the figure of merit for material characteriza- 
tion. 

Several investigators have reported the mobilities in epitax- 
ial ZnO layers grown by various techniques @, la lj] • 
The room-temperature values are relatively poor compared 
with bulk ZnO. Our recent mobilities LJl surpassed, however, 
the best bulk data ever reported and are signifi- 

cantly higher than the calculation results performed for bulk 
crystals by Rode (p. 49 of Ref. 11). Since the calculated result 
presents the material limit fill Hal , it needs to calculate the 
mobility again with a different set of parameters so as to in- 
crease its limit. The corresponding critical discussion of their 
choice will be made in this work. 

To date, there have only been a few attempts Si to model 
electron transport in doped ZnO. The present publication also 
reports carrier concentration dependence of mobility includ- 
ing an effect of the electron-plasmon interaction Q313- 

Table I gives the ZnO materials parameters that are used 
in these transport theory calculation fl4l [l5t flrill . Since the 
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current is perpendicular to the c-axis in our case Hill , the 
piezoelectric coefficient of 0.21 was used, the definition of 
which is found in p. 36 of Ref. El We have not adopted 
the relaxation time approximation for the mechanisms that in- 
volve relatively high-energy transfers, e.g., longitudinal op- 
tical phonons. The polar optical-phonon energy, Hujlo, is 
72 meV, and corresponding Debye temperature (837 K) by 
far exceeds the room temperature |16]. The pair of eo and 
too was selected from the literature as to make their differ- 
ence smaller. The effective band mass was adopted instead of 
its polaron mass. Otherwise, the resulting mobility become 
smaller. Since Rode's iterative technique takes a long time 
to reach its convergence 

Il3ll . the present computations 
are based on the variational principle method 1 14, 17]. This 
method allows the combination of all scattering mechanism 
without invoking Mattheissen's rule. 

The following electron scattering mechanisms are consid- 
ered: (1) polar optical-phonon scattering, (2) ionized-impurity 
scattering, (3) acoustic-phonon scattering through the defor- 
mation potentials, and (4) piezo-electric interactions fl4l fl5l 
Il7ll . The final expression used is as follows 0: 
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The terms eo and refer to the low- and high-frequency 
dielectric constants, respectively, m* refers to the electron 
effective mass in units of the free electron mass in rest. T 
refers to the temperature, and z% is the reduced polar opti- 
cal phonon energy defined as fojjLo/kT. Fi^iv) is the half- 
integral Fermi-Dirac integral and r\ is the reduced Fermi en- 
ergy given in thermal units. The terms -D3/23/2 an d D are 
defined and explained in Ref. 

All films were grown on insulating and lattice-matched 
ScAlMg04 (SCAM) substrates by laser molecular-beam epi- 
taxy. ZnO single-crystal and (Ga,Zn)0 ceramics targets were 
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ablated by excimer laser pulses with an oxygen flow of 1 x 
1CP 6 torr. We have grown them at about 1000 °C with intro- 
duction of semi-insulating (Mg, Zn)0 buffer layers. The de- 
tailed description of fabrication methods has been given else- 
where llSl llSffl . The films were patterned into Hall-bars and 
the contact metal electrodes were made with Au/Ti for n-type 
film, giving good ohmic contact. Carrier concentrations n s 
and Hall mobilities \i were measured at various temperatures. 

Figure 1 shows temperature dependence of mobility (^) for 
an undoped ZnO film and for a bulk crystal by closed and open 
circles. The values of // are 440 cm 2 /Vs and 5,000 cm 2 /Vs at 
300 and 100 K, respectively. Below 100 K, the distribution 
of activated carriers in this sample is too sparse to detect the 
sizable Hall electromotive force. We did not observe abnor- 
mal temperature dependence requiring a so-called two-layer 
analysis in these data 1 20] . It thus suggests the negligible con- 
tribution from a degenerate interfacial layer to our samples 
probably because ours were grown on lattice-matched sub- 
strates. 

We have derived partial mobilities by accounting for re- 
spective scattering mechanisms in the nondegenerate (un- 
doped) limit. The results are given in Fig. 1 as a function 
of temperature by dashed lines. The curve for ionized impu- 
rity scattering are calculated assuming that the material is un- 
compensated and concentration being equal to n s . The carrier 
concentration data have been given elsewhere Q|; data were 
fitted with activation energies of ^30 and 60 meV. The total 
electron mobility calculated by combining all of the partial 
scattering mechanisms is also given by a solid line. The rela- 
tive importance of the various scattering mechanisms mimics 
that found in ZnS fnll . For example, the polar optical phonon 
scattering at higher temperatures controls inherent mobility 
limit curve. At low temperatures (down to 100 K) on the other 
hand, piezoelectric phonon interaction is the dominant scatter- 
ing mechanism. Our experimental data are in reasonably good 
agreement with theory. The mobility limit at 300 K is about 
430 cm 2 /Vs. The probable discrepancy between them is be- 
cause our model calculates drift mobility, while experimental 
results correspond to Hall mobility. The theoretical Hall fac- 
tor was given in Ref.ll lL 

On the other hand, the situation is somewhat different for 
the cases of Ga-doped n-type films 12111 . Figure 2(a) shows 
300-K experimental mobilities plotted against carrier concen- 
tration. The mobilities of doped films are significantly smaller 
than those of undoped one I22L l23ll . Obviously, this ten- 
dency can be qualitatively attributed to the increased density 
of impurities. For quantitative comparison, partial mobilities 
are calculated and given in Fig. 2(a) by dashed lines. We 
have taken the effects of screening for both ionized impuri- 
ties and polar optical phonon scattering into account. Polar 
interactions reflecting the ionicity of the lattice are dominant 
in scattering mechanism, while, at heavier doping levels, ion- 
ized impurity scattering controls the inherent mobility limit 
curve I14l[l7ll . The experimental data agree well with our cal- 
culation (solid curve) except for the intermediate concentra- 
tion range. Particularly, our model could not reproduce a rela- 
tive minimum in the "// versus n s " curve, which has been ex- 
perimentally observed in this intermediate doping range. The 
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T. Makino et al., Fig. 1 (001529APL) 

FIG. 1 : (Color online). Components (dashed curves) of the total elec- 
tron mobility (solid line) at nondegenerate case versus temperature in 
ZnO. Calculations are for nondegenerate conditions. Also shown by 
closed and open circles are experimental data from an undoped film 
and a bulk crystal (cited from Ref. 10). 



situation at n s > 10 20 cm -3 could be improved probably if the 
effects of non-parabolic band structure as well as of clustering 
of charged carriers would be taken into account 0], which is 
beyond the scope of our work. 

Fischetti fl3ll and Lowney and Bennett fl4ll obtained the- 
oretical \i curve having a "dip" in the intermediate concen- 
trations by including a higher-order effect, i.e., plasmon scat- 
tering. Plasmons are the collective excitations of free elec- 
trons against the background charges. It has been shown that 
if decay due to collisional damping dominates over Landau 
damping, the electron-plasmon scattering can reduce some- 
what the mobility for electron concentrations greater than 
10 17 cm~ 3 . The coupling of the plasmon to the polar optical- 
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FIG. 2: (a) Comparison of drift mobility calculations (solid curve) 
with the Hall effect measurements for undoped and doped epitaxial 
films (filled circles). The contributions of various scattering mech- 
anisms to the total mobility are also shown by dashed curves, (b) 
Total theoretical mobility with (lower, full line) and without (up- 
per, dashed line) including the effect of electron-plasmon interaction. 
Also shown in the inset is the partial mobility. 



or emit a plasmon of energy hui p i 
fo(E') 



n s e 2 /com* is given by; 
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where upper and lower signs represent absorption and emis- 
sion, respectively. X c = [1 + * - (q c /k) 2 /2^ 1/2 ], q c is the 
plasmon cut-off wavenumber, above which plasma oscillation 
cannot be sustained. The cut-off wavenumber is on the order 
of the inverse Thomas-Fermi screening length defined else- 
where 1 14]. E' = E ± hw p i, and N p is the Bose occupation 
number of plasmons and ^ = (1 ± huj p i/E). The relaxation 
time approximation was adopted here because it is valid for 
the entire concentration range of interest according to Ref.fl3l 

The dependence of the cut-off wavenumber and plasma fre- 
quency on doping density leads to interesting phenomena for 
plasmon scattering such as a relative minimum in mobility as 
a function of concentration. Figure 2(b) shows the effect of 
electron-plasmon interaction on the mobility. The uppermost 
curve is drawn without including this effect, while the next 
one takes it into account. Notice, in particular, the "dip" is at 
n s ~ 5x 10 18 cm" 3 . The partial mobility is shown in the in- 
set. Our theory thus describes a rise due to the increase in the 
plasmon frequency at these doping levels as has been observed 
experimentally, although there is poor quantitative agreement 
between theory and experiment. We propose that more exper- 
iments using further optimized samples be done in this doping 
range to determine if the tendency on mobility agrees not only 
qualitatively but also quantitatively. 

In summary, electron transport characteristics of ZnO are 
experimentally studied and calculated by solving the Boltz- 
mann transport equation exactly for all major scattering mech- 
anisms. The polar optical phonon and ionized impurity elec- 
tron scattering mechanisms are properly screened in this treat- 
ment. Temperature dependence of the mobility is investigated 
for nondegenerate samples, and compare with our experimen- 
tal data. Reasonable agreement between experimental data 
and theory was found. Mobility limits are derived for this ma- 
terial as 430 cm 2 /V s at 300 K. Electron mobility was also cal- 
culated as a function of carrier concentration. Even if we take 
the effect of electron-plasmon interaction into consideration, 
our experimental mobilities for doped films are significantly 
lower than theoretical values especially for the concentration 
range around 10 19 cm" 3 . 
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phonon modes have been neglected because there is no suf- 
ficiently tractable theory for the hybrid modes, and because 
such coupling should not have much of an effect on the over- 
all numerical results for the mobility. We thus decide to deal 
with the pure uncoupled modes throughout. 



The momentum relaxation time t v i for an electron to absorb 



4 



[1] A. Tsukazaki, A. Ohtomo, T. Onuma, M. Ohtani, T. Makino, 

M. Sumiya, K. Ohtani, S. F. Chichibu, S. Fuke, Y. Segawa, 

et al., Nature Materials 4, 42 (2005). 
[2] A. Tsukazaki, A. Ohtomo, M. Kawasaki, C. H. Chia, 

T. Makino, Y. Segawa, T. Koida, S. Chichibu, and H. Koinuma, 

Appl. Phys. Lett. 83, 2784 (2003). 
[3] K. Miyamoto, M. Sano, H. Kato, and T. Yao, Jpn. J. Appl. Phys. 

41, L1203 (2002). 

[4] K. Elmer, J. Phys. D. 34, 3097 (2001). 

[5] D. S. Ginley and C. Bright, MRS Bull. 25, 15 (2000). 

[6] H. Kato, M. Sano, K. Miyamoto, and T. Yao, Jpn. J. Appl. Phys. 

42, 2241 (2003). 

[7] K. Iwata, P. Fons, S. Niki, A. Yamada, K. Matsubara, K. Naka- 
hara, and H. Takasu, Phys. Status Solidi (a) 180, 287 (2000). 

[8] E. M. Kaidashev, M. Lorenz, H. von Wenckstern, A. Rahm, 
H.-C. Semmelhack, K.-H. Han, G. Benndorf, C. Bundesmann, 
H. Hochmuth, and M. Grundmann, Appl. Phys. Lett. 82, 3901 
(2003). 

[9] D. C. Look, J. W. Hemsky, and J. R. Sizelove, Phys. Rev. Lett. 
82, 2552 (1999). 

[10] D. C. Look, D. C. Reynolds, J. R. Sizelove, R. L. Jones, C. W. 

Litton, G. Cantwell, and W. C. Harsch, Solid State Commun. 

105, 399 (1998). 
[11] D. L. Rode, Semiconductors and Semimetals 10, 1 (1975). 
[12] J. D. Albrecht, P. P. Ruden, S. Limpijumnong, W. R. L. Lam- 

brecht, and K. F. Brennan, J. Appl. Phys. 86, 6864 (1999). 
[13] M. V. Fischetti, Phys. Rev. B 44, 5527 (1991). 
[14] J. R. Lowney and H. S. Bennett, J. Appl. Phys. 69, 7102 (1991). 
[15] K. Seeger, Semiconductor Physics: An Introduction (Springer 

Verlag, Berlin, 2002), 8th ed. 
[16] E. Mollwo, in Semiconductors: Physics of II-VI and I- 

VII Compounds, Semimagnetic Semiconductors, edited by 



O. Madelung, M. Schulz, and H. Weiss (Springer, Berlin, 
1982), vol. 17 of Landolt-Bdrnstein New Series, p. 35. 

[17] H. E. Ruda and B. Lai, J. Appl. Phys. 68, 1714 (1990). 

[18] A. Ohtomo and A. Tsukazaki, Semicond. Sci. & Technol. 20, 
SI (2005). 

[19] T. Makino, Y. Segawa, M. Kawasaki, and H. Koinuma, Semi- 
cond. Sci. & Technol. 20, S78 (2005). 

[20] D. C. Look and R. J. Molnar, Appl. Phys. Lett. 70, 3377 (1997). 

[21] T. Makino, Y. Segawa, A. Tsukazaki, A. Ohtomo, S. Yoshida, 
and M. Kawasaki, Appl. Phys. Lett. 85, 759 (2004). 

[22] T. Makino, C. H. Chia, Y. Segawa, A. Ohtomo, K. Tamura, 
M. Kawasaki, and H. Koinuma, Phys. Status Solidi (b) 299, 
853 (2002). 

[23] T. Makino, K. Tamura, C. H. Chia, Y. Segawa, M. Kawasaki, 
A. Ohtomo, and H. Koinuma, Phys. Rev. B 65, 121201(R) 
(2002). 



TABLE I: List of material parameters. (Ref. fT^) 



Parameter 


Symbol 


Value 


Energy gap (300 K) 


E s (eV) 


3.37 


Electron effective mass 


m* (m ) 


0.24 


High-freq. dielectric constant 


too 


8.1 


Static dielectric constant 


eo 


4.0 


Optical phonon energy 


hujLO (meV) 


72 


Acoustic deformation potential 


E d (eV) 


3.8 


Longitudinal elastic constant 


Q (10 7 dyn/cm 2 ) 


2.05 


Piezoelectric coefficient 


P± 


0.21 



